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ABSTRACT 

Using Magellan/IMACS images covering a 1.2 x 1.2 sq. degree FOV with 
seeing of 0.4"-0.6", we have applied convolution techniques to analyse the light 
distribution of 364 confirmed globular cluster in the field of NGC 5128 and to 
obtain their structural parameters. Combining these parameters with existing 
Washington photometry from Harris et al. (2004), we are able to examine the 
size difference between metal-poor (blue) and metal-rich (red) globular clusters. 
For the first time, this can be addressed on a sample of confirmed clusters that 
extends to galactocentric distances about 8 times the effective radius, He//, of the 
galaxy. Within 1 R e //, red clusters are about 30% smaller on average than blue 
clusters, in agreement with the vast majority of extragalactic globular cluster 
systems studied. As the galactocentric distance increases, however, this differ- 
ence becomes negligible. Thus, our results indicate that the difference in the 
clusters' effective radii, r e , could be explained purely by projection effects, with 
red clusters being more centrally concentrated than blue ones and an intrinsic 
r e -R gc dependence, like the one observed for the Galaxy. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: individual 
(NGC 5128) — galaxies: star clusters — globular clusters: general 



1 This paper includes data gathered with the 6.5 meter Magellan Telescopes located at Las Campanas 
Observatory, Chile 
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Introduction 



Since sizes and structural parameters of globular clusters (GCs) in different GC sys- 
tems (GCSs) have first been obtained, it has become clear that some of these properties 



correlate with global p roperties of their host galaxies (see for example IJordan et al.l 12005 



Brodie Straderl 120061 ) . The existence of the so called fundamental plane relation for an 
increasing number of stu died GCSs seems to confirm that GCs populate a narrow region 
in th i s parameter space (IDjorgovskil Il995l ; iMcLaughlinl |2000| ; iMcLaughlin fc van der Marel 



20051 ; iBarmby et al. 



20071 ). However, there are puzzling trends that are still awaiting confir- 



mation and need to be addressed using larger samples of GCs. 

It is necessary to study structural parameters of GCs and GC-like objects in different 
environments before definitive statements can be made regarding their formation. Among 
the structural parameters that can be studied, the effective (or half-light) radius is of partic- 
ular importance. Model s have shown that this quantity remains fairly c onstant throughout 
the entire GC lifetime (jSpitzer &: Thuanl Il972l ; lAarseth fc Heggidll998l ). making it a good 
indicator of proto-GC sizes that are still observable today. A de cade ago, HST observation s 
unveiled a systematic size difference between red and blue GCs (IKundu fc Whitmordll998l ). 
Since then, multiple studies have found that the blue GCs are between 17% — 30% larger 



than their metal-rich counterparts in bot h spirals and early-type galaxies (IKundu et al. 
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9: 



Puzia et al 



Jordan et al 



20_05|). However, most of these studies have made use of HST observations and 



19991 : lLarsen et allboOlal Jbl: IKundu fc WhitmorelhoOll barmbv et al 



2002 



examine only the innermost regions of the galaxy or small fields in regions at galactocentric 
distances greater than the galaxy's effective radius. 



According to lLarsen fc Brodid (120031 ). the systematic size difference between red and 
blue GCs is caused merely by a projection effect. Since red (metal-rich) GCs are found to be 



more centrally concentrated than blue (meta l-poor ones) in early type galaxies (jCote et al. 



200ll ; iDirsch et al.ll2003l ; IWoodley et al.ll2005l . among others), the red GCs will appear to lie, 
on average, at a smaller galactocentric distance. The red clusters will on average be smaller 
than the blue clusters assuming that both types shares the same relation between the GC 
size and galactocentr i c dist ance. The relation r ~ iVRg C was first found in the Milky Way by 



van den Bergh et al.l (I1991I ). In this scenario, the difference between the cluster sizes should 
be most apparent at s mall galactocentric dist ance and should decrease strongly beyond 1 
galaxy effective radius (lLarsen &: Brodid 120031 ) . 



Alternatively, IJordan I (120041 ) suggests that this effect could be explained by an intrinsic 
difference between metal-rich and metal-poor GCs. Assuming half-mass radii that are in- 
dependent of metallicity, effects of mass segregation combined with a metallicity-dependent 
stellar lifetime should lead to different sizes between the blue and red clusters. The brightest 
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stars would be more massive and more centrally concentrated for the metal-rich GCs. This 
scenario should have little to no dependence on a cluster's distance from the center of its 
parent galaxy. 



In a recent study, ISpitler et al.l (120061 ) analysed the GCS of NGC 4594 (Sombrero, at a 
distance of ~ 9 Mpc) using a six-image mosaic from HST/ACS. They confirm that within the 
inner 2 arcmin (2.2 R e //), the metal-rich GCs are, on average, 17% smaller than the metal- 
poor clusters. However, the size difference becomes negligible at ~ 3 arcmin, corresponding 
to ~ 3.4 R e //, where R e // = 0.89 arcmin (IBaggett et al.lll998l ). 



To further understand the sizes of red and blue clusters, we need a homogeneous survey 
of a GCS with the ability to eliminate contaminating sources, high resolution to measure 
structural parameters, and over a large range in galactocentric distance. 



NGC 5128 is the nearest giant elliptical galaxy, at a distance of 3.8 Mpc ([McLaughlin et al 



20071 ). Its GCs are thus easily resolvable with sub-arcsecond seei ng (IHarris et al. 



this paper we present effective radius results for 337 GCs from the lWoodley et al. 



2006J). In 
2007) cat- 



alog that are confi rmed GCs by either r adial velocity measurement from variou s studies (see 



the r eferences in IWoodley et al.l 120071 ) or are resolved by HST/ACS images (IHarris et al. 



20061 ). We also present the effective radii of 27 GCs newly confirmed through radial ve- 
locity measurements using the Baade 6.5-m telescope with the instrument LDSS-2 (data in 
preparation for publication). This list represents a clean sample of confirmed clusters. All 
of these also have ellipticities less than 0.4 and effective radii less than 8 pc, both of which 
are consistent with no rmal GC prope r ties i n NGC 5128. We find that only an additional 
2.4% of GCs from the IWoodley et al.l (120071 ) catalog have effective radii greater than the 8 
pc boundary we have imposed here (to be discussed in detail in Gomez & Woodley, 2008, 
in preparation). Those few GCs are not considered here as our purpose is to establish the 
effective radius trends within the bulk of the GC population. 



2. Observations 

On the night of April 9, 2006, 25 fields were imaged with the Magellan 6.5 m tele- 
scope using the Inamori Magellan Areal Camera and Spectrograph (IMACS). In the highest 
imaging resolution, IMACS offers a FOV of 15.4 arcmin on a side, composed of a mosaic 
of 8 2Kx4KCDs with a scale of 0.111 arcsec/pixel. Our observational material will be fully 
discussed in Harris et al. 2008 (in preparation). The total field of view of our images is 
roughly 1.2 x 1.2 square degrees and the average seeing is about 0.5" across the entire field 
with individual frames ranging from 0.35" to 0.7". Images were acquired through B (on 16 
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of the 25 fields) and R (on all 25 fields) filters with both 10 second and 300 second exposures 
to avoid saturation of the brightest clusters. 



We have identified all GCs in the c atalog of IWoodlev et al.l (120071 ) on our IMACS 
framed). We have run the code ISHAPE JLarsen I [l~9~9~9l . 120011 ) individually on each GC in 
our R filtered IMACS frames, using a stellar point spread function (PSF) modelled from the 
chip in which the cluster is located. For this, typically 20-30 stars were chosen in each frame 
and measured with standard tools in ISHAPE convolves the PSF with analytical 

profiles and compares the result with the input image until a best match is achieved. 

As the analytical model, we chose King (1962) profiles, given their simplicity and because 
they are known to provide a good fit to a large family of GCs in different environments. 
Moffat functions were also tried, but they do not improve the fits except for a handful of 
large and very elliptical sources that we are not considering for the present study. They 
will be discussed as special cases in a forthcoming paper. Possible systematic ef fects in the 
sizes, arising from the choice of a particular model are discussed in lLarsen I (11999I ). However, 
the effective radius seems to be independent of the model for sources that have a similar 
extension to the stellar PSF, as in our case. For a rece nt comparison between different 



models, the reader is referred to iBarmby et al.l (120071 ) and iMcLaughlin et al.l (120071 ) . 



King I (119621 ) profiles are defined by a core-radius r c and a concentration index, which we 
define here as c = r t /r c , where r t is the tidal radius of the cluster . Usually, the concentration 



parameter is the most uncertain one to constrain (jLarsen Il200ll ). but given the high spatial 
resolution of our IMACS images, we were able to fit this along with the ellipticity, position 
angle (PA) and r c . 

The sizes quoted by ISHAPE were transformed into effective radii using the approxi- 
mation r e /r c rs 0.547c ' 486 , good to ±2% for c > 4 (jLarsen 1120011 ). The median value for 
the concentration parameter for GCs in NGC 5128 was c = 39.4 ± 10.2. Uncertainties in the 
effective radius were estimated by the standard deviation of the determined value using King 
profiles with fixed concentration indices of 15, 30, and 100. These concentration parameters 
were chosen based on typical values observed in our Galaxy as well as the concentration 
parameters fit freely with ISHAPE for the NGC 5128 data. The r e determined for any given 
cluster with varying concentration parameters does not vary more than ~ 10% for the aver- 



2 Note that the positions of two GCs have been corrected: GC0001 with a corrected right ascension of 
13' 1 25 m 1.16 s (J2000) and GC0002 with corrected declination of -43° 02 42.9" (J2000). 

3 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy Inc., under cooperative agreement with the National Science 
Foundation. 
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age GC. The concentration parameter, c, is the most uncertain of the fitted parameters. The 
extension of the GC is a secondary uncertainty. A GC at the distance of NGC 5128, with 
an effective radius of 6 pc, would span a diameter of 0.6", marginally larger than the typical 
stellar FWHM. For smaller or more compact objects, the instrisic size can be as small as 
0.1", i.e. completely blurred even with sub-arcsecond imaging. 



3. Results 



We have 69 GCs in common with the HST/ACS structural parameter study of lHarris et al 
( 120061 ). Their r e values, derived through a n isophotal analysis of the resolved clusters, and 



discussed fully in McLaughlin et al.l (120071 ) , serve as an external comparison and quality test 
for our measurements. Figure [1] shows clearly that there is good agreement in the r e values 
determined by these independent techniques. We have also examined our measured r e as a 
function of ellipticity and luminosity and also found no notable correlation. 

In an upcoming paper, (Gomez & Woodley, 2008, in preparation), we will discuss the 
structural parameters in detail as well as the new GCs discovered with the Baade 6.5m 
telescope (mentioned above) that have been used in this study. Here, we focus on the 
dependence of the GC sizes as a function of galactocentric radius, R gc for the GC subpopu- 
lations. A metallicity break has been chosen to represent the red or metal-rich ([Fe/H]> —1) 

e/H|< — 1) subp o pulations of clus t ers, f ol lowing the st u dies o f 



Larsen & Broc 


ie 


J 2003) 


Woodlev et al. 


( 


2007). 


(Harris & Harris 


2002) 



The [Fe/H] values were obtained from a C — T\ transformation 



Figures [2] and [3] show that, within one effective radius of the galaxy, the red clusters 
are significantly smaller than the blue ones by ~ 30%. As galactocentric distance increases, 
however, this difference tends to disappear and beyond a distance of 12 kpc (corresponding 
to ~ 2.3 R e //), no difference remains. We have performed a Spearman test to study the 
trend between metallicity and r e . The Spearman non-parametric rank-correlation coefficient 
ranges from -1 to 1, with -1 for a complete anti-correlation, for no correlation and 1 for 
a complete correlation. We found the coefficients to be -0.4 for K gc < 1 R e //, -0.1 for 1 
R e // < R gc < 2 R e //, 0.1 for 2R e jf < R gc < 3R e //, and 0.1 for R 9C > 3 R e //- In addition, 
the Spearman test gives the probability that the two datasets are uncorrelated as well. The 
values read 0.01, 0.19, 0.43 and 0.34, respectively, for the four radial bins presented above. 
Thus, there is a small, but virtually confirmed anti-correlation within 1 R e //. 



This size trend is even more clear in Figure [31 Within 1 R e // of the galaxy, the metal- 
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rich and metal-poor clusters have median effective radii of 1.94±0.19 and 2.98±0.20 pc, 
respectively. However, outside of 1 R e //, the median r e is identical within uncertainties. 

Figure H] shows r e as a function of projected K gc for the GCs in NGC 5128 and in 
the Milky Way. The metal-poor GCs in NGC 5128 do not follow the r e -R 5C relation- 
ship that is evident in the metal-rich GCs. The Milky Way data, on the other hand, has 
been shown to host t he r P - R gr relationship for both metallicity populations, first noted by 



van den Bergh et al.l (119911 ). using a 3-dimensional R sc . However, in projection, the metal- 



poor GCs in the Milky Way do not appear vastly different from those in NGC 5128. 



4. Discussion 



The existence of a systematic difference in the effective radii of blue and red clusters has 
been extensively studied in other galaxies, with blue clusters ty pically found t o be 17 — 30% 
larger than red ones (see Section [TJ. However, in NGC 5128, lHarris et al.l (120021 ) did not 
find any correlation b etween color and siz e for a sample of 27 GCs using HST/WFPC2. 
In a subsequent study, 1 Gomez et al.l (120061 ) found the red GCs to have larger median sizes 
compared to the blue clusters for a sample of 38 objects, with Magellan/MaglC. Both studies 
were based on small field images, involved small sample sizes, and were centered at large 
R 9C . (At the smallest R pr ,, this was mo r e than 2 times farther than the R e // of the galaxy 



light.) According to lLarsen fc Brodid (120031 ). the average sizes of red and blue clusters 
should be similar at about 1 R e f£ and beyond, if p roje ction effects a r e to a ccount for the 
size difference. The results from lHarris et al.l (120021 ) and iGomez et al.l (120061 ) are consistent 
with this scenario, bearing in mind the low number statistics of their studies. 

The only two studies thus f ar with a large enou gh sample of GCs that also extends 
beyond 1 R e //, are this work and ISpitler et al.l (120061 ). Both show that the red clusters are 
smaller than the blue within 1 R e / fj and that th e y are identical in size beyond this distance. 
In the ACS Virgo Cluster Survey, IJordan et al.l (120051 ) have studied the sizes of GCs in 67 
early-type galaxies. Their analysis reaches about 3 times the effective radius of the massive 
ellipticals studied and these GCSs dominate their sample. However, their results, favoring 
red clusters being consistently smaller than blue clusters are dominated by the inner GCs in 
these galaxies. 

Our sample, consists of only confirmed GCs which have been analyzed homogeneously, 
span a projected galactocentric distance of up to 50 kpc, i.e., 8 galaxy R e //. Thus, we are 
able to draw conclusions about the origin of the size difference with a sample that is both 
uncontaminated and much more spatially extended than in previous studies. 
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As is evident fr om Figs.[2l[3]and|l]metal-poor clusters do not show an r e -R„ c relationship. 



Jordan et al.l (120051 ) analyse this trend for metal-poor clusters in their samples and conclude 



that they are too shallow compared to the Galaxy for the projection effects to account for 
the size difference. Our results for NGC 5128 agree with this, but at the same time make it 
clear that the metal-poor subpopulation does not represent the global r e -R ffC trend. In fact, 
only metal-rich clusters show it. Therefore, projection effects can account for the observed 
size differences without the need of intrinsic formation and destruction mechanisms between 
red and blue clusters. 



5. Conclusions 

Using a contaminant-free sample of 364 GCs in NGC 5128, confirmed with radial ve- 
locity measurements or by resolved HST images, we have measured effective radii using 
ISHAPE. Our results indicate that the blue or metal-poor clusters do not show any sig- 
nificant r e -R 5C relation. However, the red or metal-rich GCs do show a steep relation in 
which red clusters within 1 R e // of the galaxy's light are 30% smaller than the blue clusters. 
Beyond this distance there is no indication for a size difference between the two metallicity 
populations. This finding in NGC 5128, not previously seen in any oth er early-type galaxy , 



supports the more tentative findings of the Sombrero galaxy's GCS ( ISpitler et al.l 120061 ) . 



Both st udies support the idea that the size differences are most likely caused by projection 



effects (ILarsen fc Brodid 120031 ) and not by intrinsic physical differences between the two 
subgroups. 
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r e (IMACS) = (1.02 ±0.04) r e (ACS) - (0.22 ±0.18)- 
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Fig. 1. — Comparison of effective radius r e , in pars ecs, for 69 GCs in NG C 5128 mea- 
sured with both our IMACS study and HST/ACS JMcLaughlin et al.l l2007h . The solid 
line is a 1:1 relationship. A least squares fit (dashed line) gives r e (IMACS) = (1.02 ± 
0.04) r e (ACS) -(0.22 ±0.18). The IMACS uncertainties correspond to the standard devia- 
tion of the effective radii in the ISHAPE fitting, using three different concentration param- 
eters c = 15, 30, 100. The agreement between these independent techniques and datasets is 
evident. 
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Fig. 2. — The effective radius, r e , as a function of [Fe/H] for (from top to bottom) R 9C < 1 
Re// of NGC 5128, 1 R e// < R gc < 2 R e// , 2i? e// < R 9C < 3R e// , and R gc > 3 R e// , where 1 
R e // = 5.1'. Least squares best fit and the number of GCs are indicated in each of the four 
panels. Only GCs with ellipticity less than 0.4 and r e < 8 pc were considered in the fit. The 
change in the slope of the fitted lines is due primarily to an increase in the median size of 
red clusters. 
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Fig. 3. — The distributions in r e of the GCs in the same galactic effective radial bins as in 
Fig. [21 for metal-poor ([Fe/H]< —1) on the left and metal-rich ([Fe/H]> —1) clusters on the 
right. Associated number of GCs and the mean r e (open circle) are labelled in each bin. 
Outside 3 ~R e ff of the galaxy, the data on the GC population suffers from incompleteness 
and spatial bias. The open circle in each distribution gives the median value for the r e and 
its formal error for each subsample. Metal-poor clusters show no significant change of r e 
with Rg C , while, on average, metal-rich clusters follow a clear trend for larger r e with R ffC . 
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Fig. 4. — The effective radius, r e in pc, as a function of projected galactocentric radius, R gc 
in kpc, for the GCs in NGC 5128 (open circles) and in the Milky Way galaxy (crosses) for 
both the metal-poor GCs ([Fe/H]< —1) on the top and meta l-rich GCs ([Fe/H]> —1) on the 
bottom. The Milky Way GC data is taken from iHarrid (11996) with the projected Ralactocen - 



tric radius defined as R gc = \/ y 1 + z 2 and R e // = 2.7 kpc (Ide Vaucouleurs k, Pence! 1 19781 ). 



Best fit curves of the form r e =c(R 5C /R e// ) Q yield (c= 3.22 ± 0.12, a = 0.05 ± 0.05) and 
(c= 2.76±0.14, a = 0.26±0.06) for the metal-poor and metal-rich GCs in NGC 5128 (dashed 
curves) along with (c= 3.02 ± 0.17, a = 0.17 ± 0.04) and (c= 3.16 ± 0.21, a = 0.36 ± 0.07) 
for the metal-poor and metal-rich GCs in the Milky Way (solid curves). 



